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Abstract 

A  two-slot  circulation  control  airfoil  was  analyzed  using  the  two-dimensional, 
compressible,  mass- averaged,  Navier-Stokes  equations.  The  implicit  Beam- Warming 
approximate  factorization  technique  was  used  to  calculate  airfoil  characteristics  for 
a  flight  Mach  number  of  0.3  and  a  Reynolds  number  near  3  million.  The  results  were 
then  compared  to  a  previous  one-slot  solution. 

An  existing  circulation  control  airfoil  was  modified  to  include  a  second  slot. 
Different  blowing  rates  were  then  applied  to  each  slot  in  various  combinations.  The 
lift  generated  for  a  given  total  blowing  momentum  for  the  two-slot  airfoil  was  nearly 
identical  to  that  for  a  single-slot  airfoil  when  the  lowest  blowing  rate  was  applied  to 
the  first  slot.  Although  the  lift  per  unit  blowing  momentum  did  not  increase  over 
the  single-slot  case,  the  maximum  lift  coefficient  was  increased  due  to  the  increased 
momentum  available  from  the  additional  slot.  Separation  angle  increased  when  a 
small  amount  of  blowing  was  applied  to  the  first  slot,  and  additional  blowing  applied 
to  the  second  slot.  The  airfoil  moment  followed  the  same  trend  as  the  single  slot, 
and  was  less  dependent  on  which  slot  the  flow  was  applied.  Due  to  the  lack  of 
experimental  data,  and  the  difficulty  in  modeling  drag  for  the  circulation  control 
airfoil,  it  is  difficult  to  compare  drag.  .  ' 
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COMPUTATIONAL  ANALYSIS  OF  A  TWO-SLOT 
CIRCULATION  CONTROL  AIRFOIL 


I.  Introduction 

Aircraft  designers  are  interested  in  high-lift  devices  for  short  takeoff  and  landing 
(STOL)  aircraft  as  well  as  rotorcraft.  One  technique  for  creating  high  lift  is  through 
airfoil  circulation  control. 

Circulation  control  incorporates  a  rounded  trailing-edge  airfoil  with  a  slot  (or 
slots)  near  the  trailing  edge.  High-pressure  air  is  blown  front  the  slot  tangential  to 
the  airfoil  surface,  creating  a  wall  jet  which  wraps  around  the  airfoil.  The  jet  ent rains 
the  airflow  from  the  upper  surface  around  the  trailing  edge  until  it  finally  separates. 
This  is  known  as  the  Coanda  effect.  The  increase  in  lift  is  roughly  proportional  to 
the  product  of  mass  flow  and  jet  velocity  (10:1).  The  stagnation  points  near  the 
front  and  rear  of  the  airfoil  shift  toward  the  center  of  the  lower  airfoil  surface.  By 
delaying  airflow  separation,  circulation  and  lift  are  increased.  Figure  1  shows  i  he 
configuration  of  a  circulation  control  airfoil.  Two  advantages  that  tin'  circulation 
control  airfoil  has  over  a  conventional  airloil  is  that  a  higher  lilt  coeffient.  C can 
be  attained  at  the  same  flight  conditions,  and  the  maximum  lift  coefficient.  Ctmax, 
is  also  increased.  This  reduces  takeoff  and  landing  speeds  for  STOL  aircraft,  and 
•ncreases  payload  for  cargo  aircraft.  Lift  can  be  varied  without  changing  the  angle 
of  attack  by  varying  the  blowing  rate.  Rotorcraft  and  X-wing  aircraft  can  benefit  by 
using  blowing  to  vary  lift  and  eliminate  the  need  for  blade  pitching  (20:1).  Because 
varying  lift  is  not  restricted  to  angle  of  attack  changes,  several  flight  controls  can  lie 
changed  by  the  blowing  rates,  such  as  climb,  pitch,  and  roll. 
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Jet  Slot 


Figure  1.  Typical  Circulation  Control  Airfoil 

The  circulation  contol  airfoil  could  replace  mechanical  flaps  but  would  require 
additional  weight  and  complexity  in  the  aircraft  wing  interior.  Aircraft  performance 
would  also  be  affected  by  energy  requirements  to  provide  airfoil  blowing  and  by  the 
negative  pitching  moments  generated  by  the  circulation  control  airfoil  exhibits. 

/./  Background. 

Many  factors  contribute  to  the  difficulty  in  the  flow-field  analysis.  The  airfoil 
lift  performance  is  very  sensitive  to  stagnation  point  locations,  therefore  the  flow 
separation  point  must  be  accurately  determined.  The  boundary  layers  of  the  flow 
after  the  wall-jets  have  different  length  scales,  and  a  free-shear  layer  forms  between 
airfoil  surface  flow  and  the  wall-jet  flow.  Compressibility  is  a  factor  even  at  low 
tree-stream  Mach  numbers  due  to  high  wall-jet  velocities. 

Several  researchers  have  applied  various  techniques  to  model  the  nonlinear  flow 
behavior  numerically  using  the  Navier-Stokes  equations.  Shrewsbury  ( 13,  15.  1-1)  and 
Pulliam  (11)  have  had  good  results  using  the  Baldwin-Lomax  turbulence  model  (3). 
with  a  curvature  correction  factor  suggested  by  Bradshaw  (5).  The  empirical  correc¬ 
tion  factored  was  adjusted  to  best  match  expermintal  data.  Williams  (ID)  used  the 
same  type  of  modeling  to  predict  airfoil  performance  for  various  turbulance  factors. 
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Harve  11  and  Franke  (8,  9)  found  experimentally  that  if  two  slots  were  used  near  the 
trailing  edge  instead  of  one,  greater  lift  can  be  obtained  using  the  same  blowing 
momentum. 

1.3  Objective 

The  purpose  of  this  study  is  lo  numerically  investigate  the  effect  of  blowing 
through  two  wall-jet  slots  compared  to  a  one  wall-jet  slot  circulation  control  airfoil. 
The  objective  is  to  evaluate  the  airfoil  performance  at  different  blowing  rates  to  de¬ 
termine  the  effect  on  lift,  drag,  moment,  and  separation.  This  is  an  extension  of  the 
work  begun  by  Williams,  and  similar  to  Shrewsbury's  work.  The  Beam- Warming 
implicit  approximate  factorization  algorithm  is  used  to  solve  the  two-dimensional, 
compressible,  mass-averaged,  Navier-Stokes  equations.  The  Baldwin-Lomax  turbu- 
lance  model  is  used,  with  the  curvature  correction  factor  by  Bradshaw  applied. 
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II.  Analysis 


The  mass-averaged,  Navier-Stokes  equations  have  been  a  very  successful  tool  in 
numerically  analyzing  the  flow-field.  The  governing  equations  for  two-dimensional, 
compressible  flow  are  described  below,  followed  by  a  description  of  the  technique  used 
to  solve  them.  The  boundary  conditions,  turbulence  model,  and  the  grid  generation 
are  then  discussed. 


3.1  Governing  Equations 

The  equations  of  motion  for  the  fluid  flow-field  are  based  on  the  conservation 
of  mass,  momentum,  and  energy.  The  variable  formulation  reflects  the  those  used 
in  the  solution  technique  developed  by  Dr.  M.R.  Visbal  (16).  In  2-D  Cartesian 
coordinates,  these  equations  appear  as  follows: 

Conservation  of  mass  (continuity) 


dp  djpu)  djpv ) 

dt  dx  dy 


=  0 


(1) 


Conservation  of  momentum 
( x-direction) 


d{pu)  djpu2)  djpuv)  Op  _  0txx  Ot yx 

dt  dx  dy  dx  dx  dy 


(y-direction) 


djpv)  djpuv)  djpv2)  dp  _  dr ^  dryy 

dt  dx  dy  dy  dx  dy 
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Conservation  of  energy 


8e  d(eu)  +  0(ev)  8(kTx )  d{kTy) 

dt  Ox  dy  Ox  dy 

d(rxxii)  8(rxy  v)  d(ryxu )  d(ryyv)  8{pu)  cj(pv) 

- - - -f-  -  -p  -  -  —  -  —  - - - 

Ox  Ox  dy  dy  dx  dy 


where: 


( 0 


Txx  =  2(/i  +  £)UX  +  A  T{UX  +  Vy) 


(<>) 


T.y  —  Tyx  —  {P  +  ~){uy  +  t?r)  l  7) 

ryy  =  2{p  +  e)vy  +  \  T{ux  +  vy)  (S) 


Other  variable  relationships  are: 

Fa  =  utxx  +  vtx  +CP(~+  -1- ) 

rr  rrt 

GA  =  UT.y  +  VTyy  +  Cp(~  +  ~) 

Stokes  hypothesis  is  applied  using: 

Ar  =  ~~{p  +  e) 

The  perfect  gas  law  is  applied  using: 


O) 

(10) 


(11) 


p  =  pRT  =  (7  -  1  )[pe  -  ^p(u2  +  v2)] 


(12) 


In  curvilinear  coordinates  (£,77),  the  Navier-Stokes  equations  can  be  combined 


as  follows: 


dq  OF  8G  OF  OG 

a7+6W  +  ^  +  ','^  +  ''‘'W  =  0 


(13) 
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where: 


/ 


F  = 


(  ,  \ 
pu 
pv 

\  Pe  ) 

pu 


pu2  +P  -  TIX 
PUV  -  TIy 

^  (pe  +  p)u  -  Fa 


( 


G  = 


pv 

pUV  —  T: 


xy 


pv2  +  p-  r. 


yy 


III) 


do) 


(ifi) 


The  transformation  between  the  physical  domain  and  the  computation;)!  do¬ 
main  is  illustrated  in  Figure  2  adapted  to  an  O-grid.  The  2-D  physical  domain 
(x.v)  and  the  computational  domain  (£,r?)  are  related  using  chain  rule  differentia¬ 
tion  (2:252). 


n  =  v(*,y) 

(1") 

d  d 

d 

(IS) 

dx  ~ 

+  Vx-^- 
or i 

d  c  0 

d 

1  Id) 

Vy  “ 

4-  r/..— 

*  (jy 

6 


Figure  2.  Physical  to  Computational  Domain  Transform 
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I  he  metrics  are  determined  bv: 


D 


d*  =  Zxdx  +  £ydy 

dr]  =  yxdx  +  rjydy 


which  in  matrix  form  are: 


1 

& 

6 

<fy  ' 

dx 

drt 

.  T]x 

Vy  _ 

dy 

uul 


dx 

d£ 

.  dy 

.  y(- 

Vr) 

.  dr]  . 

Therefore 


£x 

^x  Vy  _ 
•vhere  the  Jacobian  (J)  is: 


t  -i 


X  £  X  rj 
V(.  Vr, 


=  J 


Hv  x  r 

~’k  x< 


J  = 


d{Z,v) 

d(x,y) 


£x 

Vx  Vy 


~  £xVy  £y Vi 


or  equivalently: 


J  = 


1 


1 


d( T’V) 


xi  xv 

Vi  iJv 


(Xtl/’ 7  Xr)!J(.  ) 


(•_><)) 

i'-'l) 

(22) 

(23) 

(21) 

(25) 

(2fi) 
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The  i  evil  ting  metrics  are: 


—  {Jt}J 

7x  =  ~ViJ  Vy  = 

Sow  the  Navier-Stokes  Equation  13  can  be  recast  as: 

dq  OF  OG 

—  H - (-  -  =  0 

dt  di  dq 


F=  j(Z:F  +  tyG) 

G  =  ~j(t)xF  +  qyG) 

To  aid  in  implicit  factorization,  Equation  29  is  rewritten  as: 


dq  dE\  dE2 
dt  dq 

dVi(q,qQ  9V2  (<?,<?„)  dW^q,^) 

+  dq 
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c3urj  +  C4Vn 
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diuv  +  d2vT, 
d2Urj  +  d3Vn 
^  diuuv  +  d2(vUr,  +  uvv)  +  d3v vv  +  d^Tr,  ) 


U  and  V  denote  the  contravariant  velocities: 


U  =  (XU  +  (iyV 
V  =  T}XU  +  T)yV 
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<m<l  the  viscous  coefficients  are  6,,ct,  and  J,  : 

‘i  =  (e +  -'1(^+0 

+  -  )srsv  '  Id) 

/-» =  +  + 

»<-  'V<£  +  ><£+«;> 
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C  2 
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•> 
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2 

Hsxhy  —  “s  yhx) 
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— (m  + 
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pr  pr  )^hx  +Syhy) 

(/,  =  (m  +  +  n2y) 

,  l, 

“2  —  ~(/2  +  Z)>7r’ly  i  ■>) 

(h  =  (n  +  ^)(rll  +  ~7y) 

^4=  Cp(^:  +  +  ^y) 

J. -  Implicit  Navier-Stokes  Algorithm 

The  Navier-Stokes  solver  used  in  the  analysis  was  developed  hy  Dr.  M  R. 
V  isbal  of  the  Wright  Research  and  Deviopment  Center  (WRI)C  •)  (  ll».  17.  is).  1  he 
implicit  algorithm  uses  the  approximate  factorizaion  technique  developed  by  Beam 
and  Warming  (4)  to  solve  the  strong  conservation  Navier-Stokes  Equation  44  i  lie 
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first -order  Euler  time  differencing  scheme  in  delta  form  is  written  as: 


I  +  A  t 


dAn  d2Mn 


-  At 


di2 


I  +  At 


dBn  02Xn 

dr]  Dr,'2 


a  <r  = 


JL(El-l\-V2)n  +  -?-[E1-\\\-\Vi) 

di  dr, 


Id) 


and: 


*  n  + 1  '  n  ,  A  *  n 

q  =  q  +  Aq 


where  n  denotes  the  temporal  index  (i.e.  qn  =  q(n  At)  )  and  the  .Jacobian  matrices 


are: 


i  IS  i 


M 


0\\ 


v  _  OW2 
d<U 


Spatial  derivatives  are  discretized  using  using  second  order  central  differencing, 
while  explicit  and  implicit  damping  terms  have  been  added  to  capture  embedded 
shocks.  Because  a  steady  state  condition  is  desired,  a  variable  time  step  may  be 
used  throughout  the  grid  for  faster  convergence  (16:19-21)  .  The  mass-averaged 
Xavier-Stokes  equations  have  been  a  very  successful  tool  in  numerically  analyzing 
the  How-Held.  The  governing  equations  for  two-dimensional  compressible  How  are 
described  below,  followed  by  a  description  of  the  technique  user!  to  solve  them.  1  he 
boundary  conditions,  turbulence  model,  airfoil  configuration,  and  the  grid  generation 
are  then  discussed. 


j.  j  Boundary  Conditions 

The  conditions  on  the  physical  domain  boundaries  must  be  related  to  the 
boundary  conditions  in  the  computational  domain.  T  his  relationship  is  shown  in 
f  igure  3.  The  boundary  conditions  to  be  met  are  the  inflow  and  outflow  at  the 


O-grid  exterior,  the  matching  conditions  at  the  grid  cut  line,  and  ihe  airfoil 
with  the  wall-jet  slots. 


Figure  3.  Boundary  Conditions 

For  the  O-grid  inflow  condition,  the  air  has  the  reference  free  st  ream  value-;  at 
a  specified  angle  of  attack.  These  arc  given  by: 


u  =  cos  a 


v  -  U0 0  sin  q 


P  =  Pi c 


P  —  Poo 


For  the  outflow  conditions,  the  flow  is  subsonic  with  no  gradients  in  the  x- 
direction.  The  inflow  to  outflow  transition  occurs  near  the  maximum  and  minimum 
coordinates  of  the  O-grid  in  the  y  direction,  see  Figure  3.  The  outflow  conditions 
are  enforced  using: 


d_ 

dx 


=  o 


i  oo ) 


Along  the  O-grid  cut  line,  the  grid  is  overlapped  two  points  past  the  cut  in 
both  directions,  allowing  enforcement  of  a  periodic  boundary  condition.  On  the 
overlapped  grid  points,  the  flow  variables  are  set  equal. 


The  airfoil  has  two  types  of  boundary  conditions,  one  for  the  airfoil  surface,  and 
the  second  for  the  wall-jet  slots.  On  the  airfoil  surface,  adiabatic,  no-slip  conditions 
are  enforced  using: 

u  =  v  =  0  i  •'>()! 


d 

<)V 


=  o 


The  wall-jet  slot  is  modeled  using  assumptions  based  on  Shrewsbury's  work  (  FT  1). 
I  he  flow  from  the  wall-jet  slots  is  assumed  to  be  isentropic,  and  have  constant  total 
pressure  and  temperature  across  the  slot  exit.  The  flow  from  the  slot  is  limited  to 
less  than  or  equal  to  sonic  velocity.  The  boundary  conditions  are  enforced  using: 
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“Y  —  1  \  ^ 

T  =  T,[l  +  — - — M2 


u  =  M  (7 RT )J  cos  <p 
v  =  M  (7 RT)1  sin  <*> 


(.'/)) 
t  60 ) 

(f.l) 


Once  pressure  is  calculated  using  Equation  58,  the  Mach  number  is  determined.  For 
the  subsonic  flow.  Equations  60  and  61  are  applied  to  the  flow-field.  For  pressure 
ratios  indicating  supersonic  flow,  the  flow  is  assumed  choked  and  is  set  to  critical 
conditions. 

The  parameter  used  for  measuring  slot  blowing  is  the  blowing  momentum 
coefficient: 


9coC 


[621 


where  V'  is  the  velocity  the  jet  would  have  if  isentropically  expanded  to  the  free¬ 
st  ream  pressure.  For  two  wall-jets,  the  total  CUT  is  the  sum  of  each  slot: 


CMr  =  CM1  + 


(63) 


The  drag  coefficient  is  offset  by  the  additional  momentum  from  the  wall-jet 
by  the  following  relationship: 


Cd  —  Cdf  +  Cdp  —  CUT  ((if) 

C,if  and  Cdp  are  the  section  drag  terms  due  to  friction  and  pressure,  respectively. 
Notice  that  for  high  values  of  blowing  momentum,  C^T,  it  is  possible  to  have  negative 
values  of  drag. 
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J.j  Turbulence  Model 


The  Baldwin-Lomax  turbulence  model  (3)  was  used  this  analysis,  with  an 
empirical  correction  factor  suggested  by  Bradshaw  (5:68-71)  for  convex  surfaces. 
The  algebraic  Baldwin-Lomax  turbulence  model  is  used  because  of  its  accuracy  over 
a  wide  range  of  applications.  A  modification  for  curvature  effects  has  been  suggested 
from  theory  and  from  experiment  (5).  For  the  Baldwin-Lomax  model,  the  turbulent 
eddy  viscosity  for  the  inner  region  is: 

(Vt  ),nner  =  1,1 ''  1 

The  curvature  correction  factor  suggested  by  Bradshaw  adjusts  the  mixing  length 
(1)  by  a  factor: 

Fb  =  1  —  05  (66) 

where 


and  0  is  an  empirical  constant.  Bradshaw  suggested  that  6  <  0  <  1-1  for  a  convex 
wall  jet  and  that 

1  -  (0.05)0  <  Fb  <  1  +  (0.05)0 

limits  the  value  of  Fb ■  Williams  (19)  found  that  for  his  analysis,  a  value  of  0  =  3.2 
matched  wind  tunnel  C;  data  best.  The  value  of  0  =  3.2  was  used  in  this  analysis. 

The  boundary  layer  transition  is  based  on  the  more  stringent  of  two  conditions, 
the  first  based  on  the  Pohlhousen  method  (12)  and  the  second  based  on  adverse 
pressure  gradient.  The  most  forward  location  for  the  two  cases  is  then  used  as  the 
boundary  layer  transition  point.  This  may  produce  higher  drag  coefficients  than 
experiment  because  transition  may  not  always  occur  with  a  small  adverse  pressure 
gradient. 


2.5  Airfoil  Configurations 


The  airfoil  used  in  this  study  is  a  modified  version  of  the  103RE  circulation 
control  airfoil  designed  and  wind  tunnel  tested  by  the  David  Taylor  Xavel  Ship  lie- 
search  and  Development  Center  (DTNSRDC)  (1).  The  103RE  airfoil  is  a  modified 
ellipse  with  a  maximum  chamber  of  1%  chord  located  at  70%  chord  and  a  maximum 
thickness  of  16%  chord.  The  leading  and  trailing  edges  have  a  reduced  radius,  vary¬ 
ing  at  the  trailing  edge  from  4.6%  chord  to  2.8%.  The  airfoil  has  one  wall-jet  slot  at 
06.88%  chord  with  an  operating  slot  height  of  (0.0023)c  .  The  trading-edge  for  the 
one-slotted  103RE  is  shown  in  Figure  4.  The  research  by  Marvell  and  Era  tike  (S.  0) 


demonstrated  that  a  greater  lift  coefficient,  C/,  could  be  obtained  by  dividing  the 
blowing  momentum,  CM,  between  two  slots  near  the  Coanda  surface.  Lift  character- 
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istics  were  determined  with  the  second  slot  at  various  locations,  with  a  frce-st ream 
velocity  of  100  ft/sec.  They  found  that  the  first  slot  was  effective  in  keeping  the  flow 
attached  to  the  surface  until  reaching  the  second  slot.  The  second  slot  produced 
higher  lift  .  C/.  per  unit  blowing  momentum,  CM,  at  an  angle,  3.  of  73.5  deg  .  The 
angle  3  is  measured  clockwise  from  was  a  vertical  line  through  the  center  of  the 
circular  surface. 

The  103RE  airfoil  was  modified  to  add  a  second  wall-jet  blowing  slot  identical 
in  size  to  the  original  slot.  The  second  slot  location  was  selected  graphically  on  the 
trading-edge  at  73  deg  from  a  similarly  sized  circle,  due  to  the  variance  in  radius  of 
curvature  for  an  ellipse.  The  angle  3  for  this  point  was  70.2  deg  based  on  the  radius 
of  curvature  at  selected  point.  At  the  selected  second  slot  location,  the  existing 
airfoil  surface  between  the  first  slot  and  second  slot  was  then  translated  outward 
(0.00247)c  in  the  direction  normal  to  the  airfoil,  as  shown  in  Figure  5.  The  opening 
size  of  (0.00247)c  was  due  to  a  slot  height  of  (0.0023)c  and  a  slot  opening  lip  of 
(0.0001 7)c,  identical  to  that  of  the  original  slot.  This  shifted  airfoil  surface,  however, 
reduced  the  original  slot  height  of  the  first  slot.  The  reduction  was  determined,  and 
the  airfoil  surface  preceding  the  first  slot  was  modified  to  increase  the  first  slot  height 
to  (0.0034)c.  The  final  first  slot  height  for  the  two-slot  configuration  was  (0.0023)c 
after  the  airfoil  modifications.  The  resulting  two-slot  airfoil  trailing-edge  is  shown  in 
Figure  6. 

2.6  Grid  Generation 

The  O-grid  for  the  analysis  was  generated  using  the  interactive  grid-general  ion 
program  (INGRID)  developed  at  the  Arnold  Engineering  Development  Center  (  A  FI)( ')  (( 
Elliptical  smoothing  was  used  to  refine  the  grid  into  more  orthogonal  sections.  The 
O-grid  outer  boundary  is  a  circle  with  a  radius  of  14  airfoil  chords  from  the  airfoil  . 
Figure  7  shows  the  entire  and  Figure  8  shows  the  inner  portion.  Due  to  the  sharp 
corners  ac  the  slot  exits  ,  the  grid  is  somewhat  skewed  near  these  points,  as  seen  in 
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Figure  5.  103RE  Airfoil  Modification  -  One-Slot  to  Two-Slot 

Figures  9,  10, and  11.  The  grid  has  254  points  along  the  airfoil,  and  SO  points  normal 
to  the  airfoil.  Hyperbolic  tangent  stretching  was  used  to  concentrate  grid  points  in 
the  regions  near  the  surface  and  the  jet  slots. 

Williams’  analysis  of  the  O-grid  for  the  one-slot  configuration  found  that  a  176 
by  80  grid  was  sufficiently  dense  to  accurately  model  the  wind  tunnel  results.  The 
O-grid  for  the  two-slot  configuration  has  the  same  spacing  normal  to  the  airfoil  as 
the  one-slot  grid,  using  80  grid  points.  Along  the  airfoil  surface,  additional  points 
are  required  to  model  the  second  jet  slot.  Grid  points  were  added  between  the  slots 
and  on  the  trailing-edge  in  order  to  maintain  the  grid  density  used  in  the  one-slot 
analysis.  The  addition  of  grid  points  was  due  to  the  grid  angle  required  to  model 
the  flow  near  the  second  slot,  as  shown  in  Figure  9. 


Figure  7.  0-Grid  Overview  -  Complete  Grid 


Figure  8.  0-Grid  Overview  -  Interior  Section 
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Figure  S.  0-Grid  Detail  -  Trailing  Edge 
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III.  Results  and  Discussion 


3.1  Slot  Blowing  Conditions 

The  103RE  airfoil  has  been  extensively  wind  tunnel  tested  by  DT.YSRDC  ( 1), 
and  was  analyzed  in  Williams’  research  (19).  The  test  conditions  referenced  are  for 
Mach  number  0.3.  and  an  angle  of  attack  a  of  0.  The  wind  tunnel  test  results  for  the 
one-slot  airfoil  are  designated  as  Points  33,  35,  36,  and  38  for  various  blowing  rates. 
The  wall-jet  slot  blowing  momentum  coefficient,  Cu.  was  varied  using  a  reference  total 
pressure  in  the  slot  plenum,  and  the  measured  temperature.  The  blowing  conditions 
for  Points  33,  35,  36,  and  38  correspond  to  the  two-slot  analysis  conditions  of  A,  13, 
C.  and  D.  Point  37  was  used  to  further  investigate  lift  trends,  and  was  designated 
as  E,  even  though  the  pressure  lies  between  conditions  C  and  D. 

The  analysis  used  combinations  of  the  four  different,  blowing  rates  (A.  If.  ('. 
D)  tor  each  of  the  two  slots.  The  case  designator  for  each  two-slot  combination  is  a 
two  letter  reference,  where  the  first  letter  is  the  blowing  condition  for  the  first  slot, 
and  the  second  letter  for  the  second  slot. 

The  blowing  rate  is  varied  by  adjusting  the  total  pressure  and  temperature, 
but  C'M  will  be  affected  by  flow-field  conditions.  Williams'  analysis  had  been  done 
using  an  effective  angle  of  attach  (ae//)  and  Reynolds  number.  Re.  adjusted  for  wind 
tunnel  effects.  Both  oe//  and  Re  were  proportional  to  C/,  which  is  itself  proportional 
to  Cu.  Effective  values  were  used  in  order  to  adjust  for  wind  tunnel  effects,  allowing 
comparison  to  wind  tunnel  data  from  DTNSRDC  (1)  for  a  =  0.  Because  no  reference 
data  existed  for  the  two-slot  airfoil  conditions  being  analyzed,  linear  interpolation 
was  used  to  obtain  the  reference  (oe//)  and  Re.  The  linear  relationships  used  to 
select  the  reference  conditions  are: 

a'jj  =  -.25  -  72.3(C<ir)  f<»8) 
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Re  =  [3.12  -  1.81(Cmt  )j  106 


v.- In' re  (  \'T  is  given  by  Equation  63  ,  and  CM  values  for  Points  33.  33.  36  and  38  wore 
used.  Because  C'UJ.  is  dependent  on  the  flow-field  solution,  an  estimate  of  C  vvas 
made  using  Cu  values  from  the  one-slot  solutions.  These  estimates  of  C  varied 
from  the  final  values  by  0.033%  to  9.31%.  The  values  used  as  inputs  are  listed  in 
Table  1. 


i.i  Results 

1  he  results  of  the  analysis  at  the  various  blowing  rate  combinations  are  dis¬ 
played  in  Table  2  and  in  Figures  12  to  15. 

To  aid  in  the  result  analysis,  the  cases  were  grouped  with  those  of  constant 
ti t st  slot  blowing  conditions  (i.e.  Cases  BA.  BB,  BC.  BD).  In  Figures  12  to  15 
the  solid  lines  connect  cases  with  constant  first  slot  biowing  conditions,  allowing 
tin*  second  slot  blowing  conditions  to  change.  For  comparison,  the  dashed  lines 
in  Figures  12  to  15  connect  cases  with  constant  second  slot  blowing  conditions 
i  i.e.  Cases  AB.  BB,  CB  .  DB),  allowing  the  first  slot  blowing  conditions  to  change. 

The  single  slot  cases  are  included  as  a  reference.  The  Mach  number  profiles, 
trailing  edge  velocity  vectors,  and  Cp  vs.  x/c  for  the  two-slot,  cases  are  included  in 
the  Appendix. 

The  lift  characteristics  were  very  dependent  on  the  first  slot  blowing  moment  mu 
coefficient.  The  lift  coefficient,  Ci,  vs.  CUT  is  shown  in  Figure  12.  For  cases  without 
blowing  from  the  first  slot,  (Cases  AA,  AB),  the  flow  separated  before  reaching  the 
second  slot  and  lift  decreased  with  additional  blowing  from  the  second  slot.  Cases 
AC  and  AD  were  not  examined  because  of  flow  separation.  It  can  be  noted  that  if 
the  airfoil  has  two  slots,  blowing  may  be  required  from  the  first  slot  in  order  for  ifie 
How  to  remain  attached  until  reaching  the  second  slot. 

The  lowest  level  of  blowing  considered  for  the  first  slot  (Cases  B  A .  BB.  BC.  BD  1 
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Table  1.  Analysis  Input  Values 


Table  2.  Aralysis  Output  Values 


Case  ! 

Cr 

c, 

£ 

O 

3 

] 

J((leg) 

C  !  r 

'-I M  i  -  <'_> 

i  33  j 

0.00000 

0.032 

0.01467 

-0.01685 

7.9 

0.0000 

n/a 

i  35 

0.00956 

0.553 

0.00752 

-0.03794 

69.8 

0.00956 

u  /a,  1 

r-36 

0.018S0 

1.094 

0.00612 

mm 

0.01880 

n/a 

38 

0.03315 

-0.10150 

95.2 

0.03315 

n/a  | 

A  A 

0.00000 

-0.01438 

5.6 

0.00000 

0.00000  j 

!  AB 

0.01022 

-0.019 

0.00580 

-0.04641 

5.6 

0.00000 

0.01022  ' 

!  BA 

0.00929 

0.501 

0.00658 

-0.03516 

57.9 

0.00929 

0.00000  ■’* 
J 

BB 

0.02058 

-0.05729 

112.7 

0.01018  i  0.01040  j 

1  BC' 

0.03052 

1.609 

0.00252 

-0.07.510 

134.0 

0.01017 

0.019S6  |j 

!  BE  ! 

0.03815 

-0.08489 

163.7 

0.0 1121 

0.02691  ; 

|  BD 

0.04454 

2.353 

0.01496 

-0.09787 

173.0 

0.01153 

0.03301 

CA 

0.01852 

-0.05841 

67.1 

0.01852 

0.00000  ; 

1CB 

0.03003 

1.446 

-.00733 

-0.07710 

96.1 

0.01881 

0.01122  ! 

mm 

0.03975 

1.905 

-.00860 

-0.09810 

120.9 

0.01940 

0.02035  !' 

|  CD 

0.05263 

2.558 

0.01424 

-0.12143 

172.3 

0.01964 

0.03299 

DA 

0.03309 

0.738 

-0.01062 

-0.09165 

0.03309 

0.00000  I 

DB 

-0.02042 

-0.10658 

90.5 

0.03310 

0.01162  'j 

DC 

0.05345 

1.932 

-0.01086 

-0.13239 

109.0 

0.03316 

0. 02029 

DD  j 

0.06676 

2.353 

-0.02053 

-0.17282 

120.9 

1  0.03316  j  0.03302 

1 
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demonstrated  approximately  the  same  lift  characteristics  for  a  given  CllT  as  the  single 
slot  case.  The  lift  curve  is  extended  farther  because  a  higher  CUT  is  available  because 
of  the  two  slots.  This  may  be  an  advantage  if  slot  height  caused  choking,  as  it  does 
in  the  higher  blowing  cases. 

At  intermediate  levels  of  blowing  from  the  first  siot,  (Cases  CA.  CB,  CC.  CD), 
additional  second  slot  blowing  did  increase  C/.  However,  the  increase  in  blowing 
from  the  first  slot  was  not  as  effective  as  increasing  the  second  slot  blowing  instead. 
Consider,  for  example,  the  initial  condition  for  Case  BB.  If  additional  blowing  is 
added  to  the  first  slot  (Case  CB),  C;  increases  from  1.149  to  1.140.  with  C,,T  = 
0.03052,  shown  in  Figure  12  as  a  dashed  line  from  BB  lo  CB.  If  the  additional 
blowing  is  instead  applied  to  the  second  slot  (Case  BC),  C/  increases  from  1.1  It)  to 
1.G09.  with  CnT  =  0.03005,  shown  in  Figure  12  as  a  solid  line  from  BB  to  BC. 

At  the  highest  level  of  blowing  for  the  first  slot,  (Cases  DA,  DB,  DC,  DD),  lift 
only  increased  marginally  for  Cases  DA  and  DB,  and  actually  decreased  for  Cases 
DC  and  DD.  This  appears  to  be  due  to  flow  separation.  Additional  first  slot  blowing 
can  actually  reduce  lift  in  these  cases. 

For  comparison  to  the  one-slot  configuration,  consider  the  cases  where  no  blow¬ 
ing  is  applied  to  the  second  slot  (Cases  AA,  BA,  CA,  DA).  No' ice  in  Figure  12  that 
for  the  same  total  blowing  momentum,  CUT,  the  lift  is  less  for  these  two-slot  cases 
than  for  the  one-slot  configuration.  This  would  indicate  that  lift  and  flow  separation 
angle  for  a  circulation  control  airfoil  is  very  sensitive  to  to  the  airfoil  smoothness 
at  the  trailing-edge.  Small  protrusions  could  cause  premature  flow  separation  and 
decreased  lift. 

The  drag  characteristics  were  influenced  by  the  first  slot  blowing  momentum 
coefficient.  The  drag  coefficient,  Cj,  vs.  CMT  is  shown  in  Figure  13.  The  drag 
did  decrease  enough  to  become  negative  at  intermediate  blowing  rates,  but  drag 
increased  for  higher  blowing  momentum.  For  the  cases  in  which  drag  increased 
(Ibises  BC.  BD.  BE,  CD,  DC),  single  precision  computations  were  used.  This  may 
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have  influenced  the  drag  values  more  than  the  other  values  because  the  drag  term, 
C.i,  is  small  compared  to  the  lift  coefficient,  C/. 

The  moment  characteristics  were  dependent  on  the  total  blowing  momentum. 
The  moment  coefficient,  Cm,  vs.  CUT  is  shown  in  Figure  14.  For  all  cases,  except 
those  with  no  blowing  on  the  first  slot,  the  moment  had  a  smaller  magnitude  than 
the  one-slot  airfoil.  The  lowest  level  of  first  slot  blowing  produced  the  smallest 
magnitude  for  moment  coefficient.  The  difference  in  moment  depends  more  on  the 
total  blowing  momentum  than  on  how  the  blowing  is  divided  between  the  two  slots. 


The  flow  separation  angles  were  influenced  by  the  total  blowing  momentum 
and  how  the  blowing  was  divided  between  the  two  slots.  The  separation  angle,  .i 
vs.  C'Mr  is  shown  in  Figure  15.  For  the  same  blowing  momentum,  CMr,  the  two-slot 
configuration  had  significantly  higher  separation  angles.  The  exceptions  were  for  the 
cases  with  no  blowing  from  the  first  slot  (Cases  AA,AB).  In  these  cases,  the  flow- 
separated  immediately  after  reaching  the  first  slot.  If  the  blowing  momentum  for 
the  first  slot  increases  past  condition  B,  the  separation  angle  decreases  for  cases  with 
constant  second  slot  blowing. 
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*  *  *  *  *  1  —Slot  Airfoil 
ooqqo  2-Slot,  Slot  1  Condition  A 
QQQQP  2-Slot,  Slot  1  Condition  B 
2— Slot,  Slot  1  Condition  C 
O.fl.Q  Q.,.9  2— Slot,  Slot  1  Condition  D 


Figure  12.  Effect  of  Blowing  on  Lift  Coefficient,  (C/)  vs.  (CM) 


*****  i  —Slot 

Airfoil 

aaaaa  2  — Slot. 

Slot 

1 

Condition 

A 

Q.QQQQ  2  —  Slot, 

Slot 

1 

Condition 

B 

AAAAA  2  — Slot, 

Slot 

1 

Condition 

C 

£.0  00  0  2-Slot, 

Slot 

1 

Condition 

D 

_ 2-Slot, 

Slot 

2 

Constant 

Figure  14.  Effect  of  Blowing  on  Moment  Coefficient,  (Cm)  vs.  (CM) 


The  lift  coefficient  did  not  increase  for  the  same  total  blowing  momentum 
coefficient  for  the  two-slot  case,  which  did  occur  in  the  research  by  Ilarvell  and 
Franke  {$.  9)  with  a  different  airfoil  and  Reynolds  number.  The  lift  characteristics 
were  quite  similar  to  those  found  by  Harvell  and  Franke.  The  first  slot  blowing 
was  important  in  keeping  the  flow  attached  until  reaching  the  second  slot.  Once 
attached,  the  lift  was  increased  and  separation  delayed  primarily  by  the  second 
slot.  When  the  two-slot  case  is  compared  to  Williams’  one-slot  analysis  (19).  the 
lift  characteristics  showed  the  same  trends.  The  lift  curve  was  extended  to  higher 
maximum  lift  values,  C/mai,  because  a  higher  total  blowing  momentum  was  reached 
using  two  slots.  However,  because  of  the  different  trailing-edge  and  interference 
between  slots,  increased  first  slot  blowing  did  not  effectively  increase  lift. 

The  drag  results  were  typical  of  what  Williams  found  (19).  that  the  model 
overpredicted  the  drag.  In  cases  where  the  blowing  was  closest  to  experimental 
values,  the  drag  for  two  slots  was  less  than  a  one  slot  case  for  the  same  blowing. 

S. .{  Computer  Resources 

The  primary  computational  tool  for  the  analysis  was  the  Aeronautical  Systems 
Division  CRAY  XMP-12  through  the  sponsorship  of  Dr.  Joe  Shang  WRDC/FIMM. 
Each  run  was  approximately  2500  iterations,  with  change  in  C/  determining  con¬ 
vergence.  For  Cases  BC,  BE,  BD,  CD,  and  DC  some  computation  was  required 
at  AFIT’s  computer  facility  on  the  ELXSI-6400.  The  CRAY  computations  were  in 
double  precision,  and  averaged  3.3  x  10-6  (seconds)/[(grid  point)(iteration)j.  The 
ELXSI  computations  were  in  single  precision  with  a  performance  level  of  3.6  x  LO-3 
(seeonds)/[(grid  point)(iteration)j.  The  data  was  transferred  to  the  AFIT  computer 
systems  and  post  processed  using  the  DI-3000  software  package.  Data  transfer  and 
storage  was  difficult  to  say  the  least,  and  the  DI-3000  software  package  is  no  longer 
being  supported  at  AFIT. 


IV.  Conclusions  and  Recommendations 


l ,  /  Conclusions 

The  :Nravier-Stokes  analysis  of  a  two-slot  circulation  control  airfoil  demonstrated 
that  the  maximum  lift  coefficient,  Ctmax,  can  be  extended  over  that  of  the  one-slot 
airfoil,  but  the  lift  per  total  blowing  momentum  coefficient  did  not  increase  over 
the  one-slot  airfoil.  The  lift  coefficient, Cf,  and  angle  of  separation,  d.  show  greatest 
increases  as  additional  blowing  is  applied  to  the  second  slot,  provided  the  first  slot 
blowing  has  sufficient  blowing  to  maintain  flow  attachment  flow  between  slots.  When 
no  blowing  is  applied  to  the  first  slot,  separation  occurs  before  reaching  the  second 
slot.  Thus,  if  two  slots  exist,  some  blowing  must  be  applied  to  the  first  slot.  When 
no  blowing  is  applied  to  the  second  slot  of  a  two-slot  airfoil,  lift  and  separation  angle 
are  significantly  less  than  for  the  same  blowing  applied  to  a  one-slot  circulation 
control  airfoil.  This  indicates  that  lift  and  flow  separation  arc  very  sensitive  to  small 
protrusions  on  the  trailing-edge. 

Drag  showed  similar  characteristics  to  that  of  the  one-slot  configuration,  but  no 
conclusion  can  be  drawn  from  the  trends.  The  moment  coefficient  was  less  dependent 
on  which  slot  blowing  was  applied. 
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Recommendations 


8 


Based  on  the  analysis,  these  recommendations  follow: 


1.  Investigate  what  minimum  level  of  blowing  is  required  for  the  flow  to  remain 


,  -  -L„J  L, v.,, 

L  cl  viivu  oct  ii  O.v  iiu  uiilVi  uo  i 


*'*iicus  fi^c-ctream  conditions  affect  attach¬ 


ment. 


2.  Cse  the  NVcier-Stokes  code  to  model  existing  two-slot  airfoil  data.  The  present 
method  may  not  be  able  to  model  flow  below  Mach  numbers  of  U.25. 

3.  Analyze  various  flight  conditions,  such  as  Mach  number  or  angle  of  attack,  in 
order  to  evaluate  the  effect  on  the  airfoil  aerodynamic  characte.Ttics.  Analysis 
can  be  done  on  the  influence  of  the  effective  angle  of  attack,  and  analysis  done 
without  using  effective  values. 

1.  Investigate  the  effects  of  a  different  O-grid  generation  program,  and  if  the 
skewness  of  the  grid  near  the  slots  could  be  reduced,  improving  analysis. 

5.  Investigate  a  configuration  the  could  be  analyzed  both  experimentally  and 
computationally  for  a  two-slot,  elliptical  airfoil. 
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Appendix  A.  Cases  with  Slot  1  at  Condition  A 


Figure  16.  Case  AA  Mach  Contour  -  Overview 
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Figure  18.  Case  AA  Trailing  Edge  Velocity  Vector  Field 
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Figure  19.  Case  AA  Pressure  Coefficient  (Cp)  vs.  x/c 
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Figure  20.  Case  AB  Mach  Contour  -  Overview 
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Figure  22.  Case  AB 
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Appendix  B.  Cases  with  Slot  1  at  Condition  B 

B.  I  Case  BA  Data 


Figure  24.  Case  BA  Mach  Contour  -  Overview 


Figure  30.  Case  BB  Trailing  Edge  Velocity  Vector  Field 
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B.3  Case  BC  Data 
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Figure  34.  Case  BC  Trailing  Edge  Velocity  Vector  Field 
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Figure  35.  Case  BC  Pressure  Coefficient  (Cp)  vs.  x/c 


57 


61 


Appendix  C.  Cases  with  Slot  1  at  Condition  C 


I 


C.l  Case  CA  Data 


Figure  44.  Case  CA  Mach  Contour  -  Overview 


Figure  47.  Case  CA  Pressure  Coefficient  (Cp)  vs.  x/c 
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C.2  Case  CB  Data 
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Figure  51.  Case  CB  Pressure  Coefficient  (Cp)  vs.  x/c 
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Figure  52 


Case  CC  Mach  Contour  -  Overview 
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Figure  54.  Case  CC  Trailing  Edge  Velocity  Vector  Field 
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Figure  59.  Case  CD  Pressure  Coefficient  (Cp)  vs.  x/c 
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Figure  68.  Case  DC  Mach  Contour  -  Overview 
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